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ABSTRACT 

We demonstrate a newly developed mid-infrared (MIR) planetary nebula (PN) selec- 
tion technique. It is designed to enable efficient searches for obscured, previously un- 
known, PN candidates present in the photometric source catalogues of Galactic plane 
MIR sky surveys. Such selection is now possible via new, sensitive, high-to-medium res- 
olution, MIR satellite surveys such as those from the Spitzer Space Telescope and the 
all- sky Wide-Field Infrared Survey Explorer (WISE) satellite missions. MIR selection 
is based on how different colour-colour planes isolate zones (sometimes overlapping) 
that are predominately occupied by different astrophysical object types. These tech- 
niques depend on the reliability of the available MIR source photometry. In this pilot 
study we concentrate on MIR point source detections and show that it is dangerous 
to take the MIR GLIMPSE (Galactic Legacy Infrared Mid-Plane Survey Extraordi- 
naire) photometry from Spitzer for each candidate at face value without examining 
the actual MIR image data. About half of our selected sources are spurious detections 
due to the applied source detection algorithms being affected by complex MIR back- 
grounds and the de-blending of diffraction spikes around bright MIR point sources into 
point sources themselves. Nevertheless, once this additional visual diagnostic checking 
is performed, valuable MIR selected PN candidates are uncovered. Four turned out 
to have faint, compact, optical counterparts in our Ha survey data missed in previ- 
ous optical searches. We confirm all of these as true PNe via our follow-up optical 
spectroscopy. This lends weight to the veracity of our MIR technique. It demonstrates 
sufficient robustness that high-confidence samples of new Galactic PN candidates can 
be extracted from these MIR surveys without confirmatory optical spectroscopy and 
imaging. This is problematic or impossible when the extinction is large. 

Key words: (ISM:) planetary nebulae: general - (ISM:) - infrared radiation extinction 
- radio continuum - Hll regions 



1 INTRODUCTION 

We present an investigation into the potential of mid- 
infrared (MIR) survey data from the Spitzer and WISE 



space satellite missions as a tool to uncover planetary neb- 
ula (PN) candidates that would be hard or impossible to 
locate optically. The motivation is to develop MIR PN can- 
didate selection techniques that can be used to uncover the 
significant numbers of Galactic PNe which are believed to 
be hidden behind extensive curtains of dust. 
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PNe are important astro-physical objects and key win- 
dows into late stage stellar evolution. They play a major role 
in Galactic chemical evolution (Dopita et al. 1997; Karakas 
et al. 2009), return significant enriched mass to the ISM 
(Iben 1995) and are powerful kinematic tracers due to their 
strong emission lines (eg. Durand, Acker & Zijlstra, 1998). 
Over the last decade Galactic PNe discoveries have entered 
a golden age due to the advent of narrow-band Galactic 
plane surveys of high-sensitivity and resolution. This has 
been coupled to complementary, multi-wavelength surveys 
across near infrared (NIR), MIR and radio regimes in par- 
ticular from both ground and space-based telescopes. These 
have provided powerful diagnostic and discovery capabili- 
ties (e.g. Cohen et al. 2007, 2011, hereafter Papers 1 and 2; 
Phillips and Ramos-Larios 2008; Ramos-Larios et al. 2009; 
Miszalski et al. 2011; Anderson et al. 2012). 

The total number of known Galactic PNe is cur- 
rently ~3000, double what it was a decade ago. This is 
largely due to the ^1200 PNe found by the two Mac- 
quarie/ A AO/Strasbourg Ha PNe surveys (MASH: Parker 
et al. 2006; Miszalski et al. 2008). MASH PNe were un- 
covered via scrutiny of the sensitive, arcsecond resolution 
SuperCOSMOS AAO/UKST Ha survey of the Southern 
Galactic plane (SHS: Parker et al. 2005). These are now 
being supplemented by equivalent discoveries in the North- 
ern Galactic plane (e.g. Mampaso et al. 2006; Viironen et 
al. 2009a, b; Sabin et al. 2010) arising from careful searches 
of the Isaac Newton Telescope Photometric Ha Survey data 
(IPHAS: Drew et al. 2005). However, these combined num- 
bers still fall a factor of two short of even the most conserva- 
tive Galactic PN population estimates (Jacoby et al. 2010) 
where population synthesis yields 6,600-46,000 PNe depend- 
ing on whether the binary hypothesis for PN formation is 
required (e.g. De Marco 2009). Significant numbers of PNe 
are faint and highly evolved as shown to exist in the local 
volume sample of Frew & Parker (2006) and Frew (2008). 
They rapidly become undetectable at distances greater than 
a few kpc. Such objects currently remain beyond detect abil- 
ity. 

However, there are also serious problems with obtaining 
truly representative samples of PNe across the galaxy due to 
variable extinction. It is clear that a significant population 
of Galactic PNe is lurking behind the extensive clouds of gas 
and dust that obscure large regions in the optical regime. It 
is the extension of previous PN discovery techniques away 
from the optically dominant [OIII] PN emission line in un- 
reddened spectra to the longer wavelength Ha emission line 
(that can peer at least partially through the dust), that has 
led to the major, recent discoveries. Extension of PN iden- 
tification techniques to longer, more favourable wavelengths 
would clearly be advantageous. 

For this pilot study six MIR colour-colour selection cri- 
teria were simultaneously applied to the 49 million entries 
in the GLIMPSE-I point source archive. These criteria are 
based on sources within three standard errors of the median 
values of the six unique MIR colours of the 136 previously 
known PNe that fall within the GLIMPSE-I footprint. These 
are assumed representative of the overall Galactic PN pop- 
ulation as given in Paper 1. Only 70 candidate sources were 
returned. About half turned out to be spurious once the MIR 
image data were examined. Despite significant extinction, 
four of the remaining sources had faint optical detections 



in the available Ha survey images and are the basis of the 
optical spectroscopic follow-up. Most Galactic PNe are well 
resolved and will not be found in searches of the GLIMPSE- 
I point-source archive which also has very restricted Galac- 
tic latitude coverage. These factors substantially reduce the 
number of obscured PN candidates returned here. 

Section 2 gives the background importance and context 
of this MIR PN study. Section 3 briefly describes current 
knowledge of PN characteristics at non-optical wavelengths 
and the importance of eliminating mimics. Section 4 de- 
scribes the candidates' MIR selection. Section 5 gives our 
new study into MIR false-colour imagery as a powerful di- 
agnostic tool. Section 6 presents the spectroscopic follow-up 
of the four optical counterparts to our MIR selected PN 
candidates. They are all confirmed as PNe. An additional, 
serendipitous PN found adjacent to one of our optically de- 
tected MIR sources is also confirmed. Some basic character- 
istics of the new PNe are also presented. In section 7, we 
provide some conclusions and suggestions for future work. 



2 PN DISCOVERIES AT NON OPTICAL 
WAVELENGTHS 

Jacoby & Van de Steene (2004) undertook early work with 
an on-band, off-band [SIII] 9532 A emission line survey in 
the Galactic bulge as this is a prominent PN line in the 
far red. They found 94 candidate PNe though many still 
require confirmation. More generally, significant PN can- 
didates have been selected via their Infrared Astronomical 
Satellite (IRAS) colours but confirmatory success rates have 
been low, compounded by the large IRAS error ellipse (e.g. 
Suarez et al. 2006, Ramos-Larios et al. 2009). This is an inef- 
ficient technique not considered further. More recently, MIR 
space-telescope images from Spitzer (Werner et al. 2004) and 
now WISE (Wright et al. 2010) allow the detection of very 
reddened PNe that may be invisible optically (eg. Cohen et 
al. 2005, Kwok et al. 2008; Phillips & Ramos-Larios 2008; 
Zhang & Kwok 2009; Zhang, Chih-Hao & Kwok. 2012). 
Carey et al. (2009) and Mizuno et al. (2010) noted 416 com- 
pact but resolved (< 1 arcmin) ring, shell and disk-shaped 
sources in the Galactic plane in 24/xm Spitzer MIPSGAL im- 
ages (MIPSGAL is an extensive infrared survey of the Galac- 
tic plane using the Spitzer Multiband Imaging Photometer 
(MIPS) instrument, see Rieke et al. 2004). Based on expe- 
rience we think many will be strongly reddened PNe with 
only a minority being circumstellar nebulae around massive 
stars (Wachter et al. 2010; Gvaramadze et al. 2010). 

PNe can be strong NIR and MIR emitting objects. This 
is because of their PAH emission, fine structure lines like 
[OIV] at 25.89/xm (e.g. Chu 2012), thermal dust emission 
within the nebulae and from circumnuclear disks and H 2 
molecular lines (e.g. the UKIRT Wide-field Infrared Survey 
for if 2 , (UWISH) Froebrich et al. 2011). Papers 1 and 2 
analysed optically detected known PNe and PN candidates 
in the Spitzer GLIMPSE-I survey (Benjamin et al. 2003) to 
develop robust, multi- wavelength classification and diagnos- 
tic tools that provide purer PN samples in heavily obscured 
regions. GLIMPSE-I is the 2-degree wide Spitzer mid-plane 
survey. The goal is to recognise quality PN candidates solely 
using MIR and radio characteristics, enabling the search for 
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hidden PNe when heavy extinction prevents optical spectra 
and images. 

2.1 Eliminating non-PN contaminants 

Non-PN contaminants have badly undermined the integrity 
of pre-MASH PN catalogues. Objects with extended emis- 
sion can often masquerade as PNe. These include compact 
HII regions (Cohen et al. 2011), Stromgren zones (Madsen 
et al. 2006; Frew et al. 2010), ejecta shells around Wolf- 
Rayet and other massive stars (e.g. Marston 1997; Chu 2003; 
Stock & Barlow 2010), supernova remnants (e.g. Stupar et 
al. 2007, 2011), symbiotic systems (e.g. Corradi 1995; Cor- 
radi et al. 2008, 2010), Herbig-Haro objects and their kin 
(Canto 1981), as well as nova shells and bow-shock nebulae 
(Frew & Parker 2010). Identification is further complicated 
by the variety of morphologies, ionization properties and 
surface-brightness distributions exhibited by the PN family 
itself. We have tested and developed criteria to more effec- 
tively eliminate contaminants using new multi-wavelength 
surveys combined with emission line ratios from follow-up 
spectroscopy. This has enabled clear discrimination tools to 
be developed (see Frew & Parker 2010). 

In Paper 2 we applied these criteria to our MIR samples 
of known optically detected Galactic PNe seen in GLIMPSE- 
I and that overlap with the SHS (i.e. \b\ ^1 degree and from 
210 through 360 to 40 degrees in Galactic longitude). This 
showed that 45% of previously known pre-MASH PNe are 
HII region contaminants. The MASH contaminant fraction 
was only 5% in the same zone as these discrimination tech- 
niques had already been applied to MASH. Furthermore, 
external filaments, structures and/or amorphous halos seen 
for MIR sources generally indicates an HII region. This is an 
important MIR diagnostic for discriminating resolved MIR 
HII regions from objects such as PNe as both object types 
can look similar optically. In a similar vein, Anderson et al. 
(2012) use IR data from Herschel (Hi-Gal), WISE, MIPS- 
GAL and GLIMPSE to independently establish IR selection 
criteria to distinguish between HII regions and PNe. 



3 SELECTION OF PN CANDIDATES FROM 
THE GLIMPSE-I ARCHIVE 

We want robust MIR selection criteria to identify quality PN 
candidates in highly obscured Galactic plane regions where 
the current PN number density is unsurprisingly low. The 
zones of avoidance in optically identified PN samples are 
due to extinction, particularly in the bulge (e.g. see Fig. 6 
of Miszalski et al. 2008) but the PNe must be there. 

Paper 1 offered three MIR colour-colour planes designed 
to isolate the PN domain from diffuse and compact /ultra- 
compact HII regions which are the dominant PN catalogue 
contaminants. PN median colours were derived from the 136 
known PNe in the GLIMPSE-I survey. We use our integrated 
photometry for each source (resolved or compact) from the 
Spitzer Infrared Array Camera instrument (IRAC; Fazio et 
al. 2004) at 3.6, 4.5, 5.8 and 8.0 /mi. These bands provide 
six colour indices examined for clear trends and locii (e.g. 
see Fig. 9, 11 and 12 in Paper 2). Although most of the 136 
known PNe are extended, for this pilot study we concentrate 
on applying our MIR PN selection criteria to entries in the 



Table 1. The 6 GLIMPSE-I MIR colour-colour selection crite- 
ria applied to the point-source archive based on median colours 
±3 sems of the 136 known PNe in the GLIMPSE-I footprint. 



IRAC 


median 


lower 


upper 


colour 


colour 


bound 


bound 


[3.6] - [4.5] 


0.81 


0.57 


1.05 


[3.6] - [5.8] 


1.73 


1.43 


2.03 


[3.6] - [8.0] 


3.70 


3.37 


4.03 


[4.5] - [5.8] 


0.86 


0.56 


1.16 


[4.5] - [8.0] 


2.56 


2.23 


2.89 


[5.8] - [8.0] 


1.86 


1.65 


2.07 



GLIMPSE-I archive which includes only point-like sources 
resolved by Spitzer/IRAC. Prospects for uncovering new re- 
solved MIR PN candidates (e.g. as found by Mizuno et al. 
2010) will be the basis of a separate paper. 

Using our median IRAC PN colours from paper 1 and 
adopting the standard errors of the medians (sem) we used 
the NASA Infrared Science Archive (IRSA) search tool to 
query all 49 million point sources in the GLIMPSE-I spring 
2007 archive. This is deeper, but with slightly larger pho- 
tometric uncertainties, than the 31 million sources in the 
GLIMPSE-I catalogue. Objects whose IRAC colour indices 
simultaneously satisfy all six colour-colour selection criteria 
and therefore fall within the median ±3 sem boxes defined in 
Table 4 in Paper 1 were selected. This applies the maximum 
rigour to our search for new PNe and the greatest rejection 
of HII regions. The resulting MIR colour selection returned 
70 PN candidates. 

Such a search does not return all 136 known PNe in 
Paper 2 nor all possible unknown PNe in the GLIMPSE-I 
footprint. There are four main reasons. Firstly, known PNe 
are usually extended and will not appear in the GLIMPSE-I 
point-source archive. Secondly, due to heavy extinction close 
to the Galactic plane, many MIR point-source candidates 
have no optical counterpart so are unlikely to have been 
previously identified. Thirdly, most known PNe (^80%) 
fall outside the median ± 3 sem MIR colour-colour boxes. 
Fourthly, not every PN has a signature or reliable photo- 
metric measurement in each IRAC band. Also only 60% of 
known PNe can be separated from their local background 
by false colour or colour-index (Paper 2) while 40% lack 
the data to define false colour or have only poor quality 
photometry. Of the 17 known PNe that fall within the me- 
dian ± 3 sem of the [3.6]-[4.5] versus [5.8]-[8.0] MIR colour- 
colour box shown in Fig|3] 69% are compact with major-axes 
<10 arcseconds. We do expect to recover the most com- 
pact/barely resolved known PNe that fulfil our criteria. 

3.1 Matches with known PNe and uncovered 
positional errors 

We searched for counterparts within 30 arcseconds of the 
70 MIR-selected PN candidates using SIMBAD (Wenger et 
al. 2000). In this way resolved sources (like PNe) that may 
have compact MIR cores not coincident with the centre of 
the nebula would be selected. Two candidates are associated 
with known PNe: Hen 2-84 (ESO 95-9; PN G300.4-00.9), a 
relatively compact bipolar PN ~30 arcseconds across and 
K 3-42 (PNG 056.4-00.9) a very compact, point-like PN 
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^3.5 arcseconds across. The bipolar PN is seen as a MIR 
point-source because of a bright, coincident MIR source per- 
haps due to a compact, heated (circumstellar) dust torus 
close to the central star or even a late-type binary compan- 
ion. It is also seen in the SuperCOSMOS Sky-Survey (SSS: 
Hambly et al. 2001) I-band data. We have assumed through- 
out that SIMBAD offers a reasonably complete inventory of 
likely catalogued cross-identifications for our candidates. 

Kerber et al. (2003) provide incorrect J2000 co- 
ordinates for K 3-42 of 19h39m36.0s +20° 19' 07" . The actual 
Ha source in the astrometrically calibrated IPHAS and SHS 
surveys is at 19h39m35.8s +20°19' 02", ~7 arcseconds away 
(though still matched with the true GLIMPSE-I source by 
the larger initial SIMBAD search radius). Given the source's 
compact nature and high stellar density this offset is signif- 
icant. The literature position gives an incorrect match to 
the GLIMPSE-I source SSTGLMC G056.4034-00.9035 for 
K 3-42. The true GLIMPSE-I source and the one with the 
appropriate MIR colours is SSTGLMA G056.4016-00.9033 
(< 1 arcsecond from our corrected position). 

Pe 2-8 (ESO 177-3) is a known, compact PN in the list 
of 136 but was not returned from the GLIMPSE-I search. 
SIMBAD returns the associated GLIMPSE-I source SST- 
GLMC G322.4689-00.1778 but it only passes three of the six 
the MIR point-source colour-colour criteria applied. None of 
the more abundant, well resolved, and/or nearby PNe regis- 
ter as a point-source in the GLIMPSE-I archive unless they 
contain an associated MIR star (as for Hen 2-84). Neverthe- 
less, these eventual matches between the MIR sources and 
their true optical counterparts (once accurate co-ordinates 
are used) show that the method has successfully retrieved 
two out of the three compact known PNe that were likely 
to be in the GLIMPSE-I point source catalogue. 



3.2 The MIR selected PN candidates 

Among the 70 returned MIR candidates 29 (42%) have no 
current catalogued identification within 30 arcseconds in 
SIMBAD. Reducing the search radius to 2 arcseconds, a rea- 
sonable choice given the point-source nature of GLIMPSE- 
I, the number of uncharacterised sources rises to 45 (64%). 
Many of these unidentified sources may be PNe. The 25 
(36%) other sources with a SIMBAD entry provide matches 
with two known PNe (as above), a YSO and 22 associated 
with suspected young stellar sources (designated as Y*? en- 
tries in SIMBAD but referred to as YSO? hereafter). How- 
ever, the status of these suspected YSOs is debatable. They 
form part of the 11,000 likely YSOs in the flux-limited cen- 
sus of 18,949 point sources in the Galactic mid-plane selected 
from GLIMPSE-I and GLIMPSE-II from their intrinsically 
red MIR colours (Robitaille et al. 2008). However, in their 
theoretical grid of 200,000 YSO spectral energy distributions 
(SEDs), Robitaille et al. (2006) offered NIR and MIR colour- 
colour plots with the locations of these SEDs. In particular 
their Fig. 18 shows the [3.6]- [4.5] indices for the full grid. 
The range of YSOs spans several magnitudes and almost all 
YSOs exceed the [3.6]- [4.5] colour of PNe. We now investi- 
gate the veracity of these identifications given that we have 
found that they appear to have the MIR properties of PNe 
according to our selection criteria. Despite any associations 
of our 70 MIR selected sources with catalogued sources of 



mostly unproven type, a fresh evaluation of their true nature 
was also needed. 



3.3 MIR PN candidates with optical counterparts 

The MIR PN candidate selection is validated if some are 
shown to be PN. As many MIR candidates lie in obscured 
regions it was important to search for optical counterparts in 
the Ha survey images that would permit confirmatory op- 
tical spectroscopy. Consequently, false colour RGB optical 
images (2x2 arcmin) of each of the 70 sources were created 
using the on-line SHS data with Ha as red, matching short- 
red (SR) as green and the broad-band SSS B j image as blue. 
Quotient images from dividing the Ha image by its match- 
ing SR counterpart were also created. These images high- 
light Ha emitters. Four objects with detectable Ha emis- 
sion at the MIR position were identified, indicating an opti- 
cal counterpart. False-colour 1.5x1.5 arcmin images of these 
four optically detected candidates are given in Fig. [I] The 
left panels comprise the optical false-colour images. Note 
the faint ness of the sources in Ha. The middle panels are 
the Ha/SR quotient images which clearly reveal the com- 
pact emitting sources. The right panels are the IRAC432 
colour combination with RGB = [8.0], [5.8] and [4.5] /mi. 
The similar orange-red false colours for all PN candidates in 
the IRAC432 images are evident. 

3.4 Serendipitous PN 'discovery' next to optical 
MIR PN candidate GLIPN1557-5430 

A compact Ha source was noticed in the 2x2 arcmin 
RGB Ha/SR/Bj finding chart 61 arcseconds away from 
GLIPN1557-5430, one of the MIR PN candidates with an 
optical counterpart. Very few Galactic PNe have a compan- 
ion within 1 arcminute (Parker et al. 2006). This source was 
previously identified as a possible PN based on its far-IR 
IRAS colours. It is listed as object 10 in Table 2 of Phillips 
& Ramos-Larios (2008) and listed as a post-AGB/PN candi- 
date in Table 1 of Ramos-Larios et al. (2009) with an IRAS 
identification and designated PM 1-104. The SIMBAD co- 
ordinates are from these references but it has never been 
confirmed as a PN until now (see below). PM 1-104 is not 
among the 70 selected MIR sources as it satisfies only 3 
of the 6 colour-colour selection criteria. Many known PNe 
fall outside the ±3 standard deviations from the median 
colour-colour selection criteria though those that fall inside 
are more likely to be true PNe as this is the purest MIR 
colour-colour space for known PNe. The source is compact, 
of relatively high surface brightness optically and has very 
similar MIR false colours to true PNe. Furthermore, like for 
K 3-42, the published position of Ramos-Larios et al. (2009) 
is erroneous by ^8 arcseconds. 

Our optical and MIR colour montages reveal the source 
to be at 15h57m21.0s -54°30 / 46 // (J2000) and not at the 
15h57m20.4s -54°30 / 40 // position reported in the literature. 
We included this object in our list for spectroscopic observa- 
tion. The MIR IRAC colours for the 4 new optically coun- 
terparts together with PM 1-104 and the two known PN in 
the MIR sample are given in Table [2] None of these optically 
detected sources (except PM 1-104) is previously known or 
recorded in MASH because they are very faint, compact and 
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Figure 1. Left panels are 1.5 x 1.5 arcminute colour composite SuperCOSMOS images comprising the Ha image (red), matching 
broad-band SR image (green) and Bj-band image (blue) of confirmed PNe GLIPN1833-1133, GLIPN1530-5557, GLIPN1557-5430 and 
PN GLIPN1642-4453. Each image is centred on the faint PN as seen in Ha. Middle panels are quotient images of Ha divided by the 
matching broad-band short red (SR) exposure that reveals the compact emitting sources. Right panel: same area from the IRAC data 
with RGB=8.0, 5.8 and 4.5/im - i.e. the IRAC432 false-colour combination. All PN candidates have similar IRAC false-colours. 



extremely hard to find. The remaining sources do not appear 
to have any detectable Ha emission. 

We have uncovered two significant published positional 
errors for PNe/possible PNe in this small sample. We are 
completely revising the published co-ordinates for all known 
and MASH PNe based on new multi-wavelength imaging 
with decent astrometry to provide the definitive list of ac- 
curate PN positions (Parker et al., in preparation). We have 



shown that the situation reported above is common with 
literature positions for over 100 PNe out by >10 arcseconds 
compared to our own accurate values based on the latest 
surveys and accurate astrometric solutions. 
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Table 2. IRAC colours for the four new MIR selected PN candidates with optical counterparts. The serendipitously uncovered source 
PM 1-104 (now also a confirmed PN) and two known PNe in the MIR sample are also included. All IRAC false-colours for these sources 
including those in Fig. ^appear similar (Orange-red). 



New PN ID 


[3.6]-[4.5] 


[3.6]-[5.8] 


[3.6]-[8.0] 


[4.5]-[5.8] 


[4.5]-[8.0] 


[5.8]-[8.0] 


GLIPN1530-5557 


0.710 


1.718 


3.409 


1.008 


2.699 


1.691 


GLIPN1557-5430 


1.038 


1.822 


3.776 


0.784 


2.738 


1.954 


GLIPN1642-4453 


0.753 


1.741 


3.509 


0.988 


2.756 


1.768 


GLIPN1823-1133 


0.891 


1.643 


3.699 


0.752 


2.808 


2.056 


PM 1-104 


0.454 


1.537 


2.948 


1.083 


2.494 


1.411 


Hen 2-84 


1.021 


1.786 


3.676 


0.765 


2.655 


1.890 


K 3-42 


0.857 


1.644 


3.631 


0.787 


2.774 


1.987 



4 THE DIAGNOSTIC POWER OF MIR FALSE 
COLOUR IMAGERY 

In Paper 1 we showed that PNe have only three basic colours 
in our combined IRAC band false colour images: red, orange 
and violet. This is due to the combination of PN narrow 
atomic and molecular lines, PAH bands of modest width and 
sometimes broad dust emission that comes through these 
bands, giving PNe their distinct hues and providing a valu- 
able visual diagnostic. Similarly, 2MASS J,H,Ks false colours 
of PNe are violet or pink/purple when they are detected. 
The advent of the Vista Variables in the Via Lactea (VVV) 
NIR survey in the Y,Z,J,H,Ks photometric bands (Saito et 
al. 2012) enables us to search for these MIR selected PN 
candidates in the NIR afresh. VVV has the same J,H,Ks 
pass-bands as 2MASS, but depth and resolution are far su- 
perior and extend ~ 4 mag deeper than 2MASS; Saito et al. 
2012). Miszalski et al. (2011) used the equivalent Vista data 
for the Large Magellanic Cloud (the VMC survey, Cioni et 
al. 2010) to show that many PNe can be seen in these high- 
quality NIR data. They showed that by constructing SEDs 
across a broad wavelength range from extant wide-field sky 
surveys it is possible to discriminate among different astro- 
physical sources including PNe. 

As an example we present in Fig. [2] mult i- wavelength 
1.5x1.5 arcminute false-colour montages of PM 1-104 and 
GLIPN1 577-5430 which, given their angular proximity, both 
appear in these extractions which are centred mid-way be- 
tween the two sources. This figure demonstrates the power of 
mult i- wavelength, false-colour imagery applied to PNe. The 
nature of the colour images are given at the bottom right 
of each image. The literature position for PM 1-104 is seen 
offset from the brighter source itself, revealing the noted 
positional error. In the VVV J,H,Ks image both PNe are 
pink sources, similar to resolved PNe in 2MASS. In IRAC 
they both appear as almost identical yellow-orange sources 
in complete contrast to every other source. In the poorer res- 
olution WISE image-band combinations they are also easily 
detected, appearing distinctly red (when combining bands 
321) or yellow (when combining bands 432). 

In Fig. [3] we present the IRAC [3.6]-[4.5] versus [5.8]- 
[8.0] colour-colour plot (equivalent to Fig. 15 in Paper 1 and 
Fig. 9 in Paper 2) with the four newly confirmed PNe (refer 
to section on spectroscopy) plotted as red circles within the 
MIR PN selection box. This particular colour-colour plane 
has the greatest discriminatory power with the ±3 sem PN 
box being furthest from the majority of other astrophysical 
source types. There is no possible confusion of these MIR 



selected point sources with the box occupied by diffuse HII 
regions which slightly overlaps the PN box. The newly con- 
firmed PN PM 1-104 is plotted as a blue symbol. In Table|3] 
we also present the 2MASS NIR and GLIMPSE-I MIR pho- 
tometric data for all of these sources. 



4.1 MIR photometry and environment: 

identification of high quality PN candidates 

All the SHS and IRAC432 false-colour images for the 70 MIR 
sources were carefully examined during the search for opti- 
cal counterparts. It became clear that not all MIR sources 
satisfying the 6 colour selection criteria were as similar in 
their false colour imagery as expected. This is due to their 
complex MIR emission environments. Many sources were 
matched in SIMBAD with a 2 arcsecond search radius to sus- 
pected YSOs from Robitaille et al. (2008). For many of these 
supposed MIR point-sources the false colours, surrounding 
MIR diffuse environment or conversely match to stars in 
the NIR 2MASS and SSS optical images strongly mitigate 
against PN identification. Where sources are embedded in a 
MIR diffuse emission region (with prominence varying across 
IRAC bands) the photometry and MIR colours have been 
adversely affected by the background. Hence, as part of our 
PN candidate selection we adopted the additional selection 
criteria: 1) There should be no obvious star co-incident with 
the MIR position in the optical and 2MASS; 2) The MIR 
environment should be as clean and unstructured as possi- 
ble; 3) The MIR source should be a compact source with 
the standard IRAC432 red-orange colour normally seen for 
all known PNe (and indeed the four newly confirmed PNe). 

With these criteria the PN candidate landscape changes 
considerably and we select only 33 of the 70 MIR sources as 
PN candidates (of which two are known PNe and four are 
the new PNe presented here). We reject the other 37 sources 
as contaminants. These 33 remaining sources have similar 
false colours to known PNe with 14 having no SIMBAD 
entry. Suspected YSO sources in SIMBAD from Robitaille 
et al. (2008) are associated with 13 of these sources but we 
contend they are in fact excellent PN candidates. 

This clearly demonstrates the danger of both trusting 
the SIMBAD identifications and over-interpreting the IRAC 
point-source photometry when used in isolation. The value 
of eye-ball examination of false colour MIR and NIR im- 
ages of each source to compare with the optical is clear. 
In Fig. [4] we present MIR 2MASS J,H,Ks, optical SHS and 
IRAC432 false colour 1.5x1.5 arcmin images of a selection 
of three of the 37 rejected PN candidates as examples of the 
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Figure 2. Multi-wavelength 1.5x1.5 arcminute false-colour montages of newly spectroscopically confirmed PN PM 1-104 (GLIPN1557- 
5431) and GLIPN1557-5430. Image type is indicated at the bottom right of each panel. The published position of PM 1-104 is identified 
by tick marks towards the lower left-hand corner close to the actual source (also marked) revealing the positional error. Both PNe are 
also visible in the two lower resolution WISE false colour images. See text for further details. 



Table 3. 2MASS and IRAC magnitudes (with la errors) for the MIR selected sources, the serendipitous recovered PN PM 1-104 and 
the two known PNe that satisfy the MIR selection criteria. 



New PN ID 


J 


H 


Ks 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


GLIPN1530-5557 


12.92±0.04 


11.59±0.04 


10.95±0.04 


10.00±0.04 


9.29±0.05 


8.28±0.04 


6.59±0.03 


GLIPN1557-5430 




14.62±0.10 


13.28±0.06 


11.77±0.04 


10.73±0.05 


9.95±0.04 


7.99±0.03 


GLIPN1642-4453 






13.74±0.07 


12.21±0.10 


11.46±0.09 


10.47±0.09 


8.70±0.10 


GLIPN1823-1133 








12.74±0.08 


11.85±0.10 


11.094±0.09 


9.04±0.05 


PM 1-104 


14.17±0.06 


13.27±0.08 


12.01±0.04 


9.20±0.12 


8.75±0.10 


7.67±0.07 


6.26±0.13 


Hen 2-84 








12.25±0.07 


11.23±0.08 


10.47±0.08 


8.58±0.07 


PN K 3-42 


14.20 ± 0.04 


13.71±0.05 


12.70±0.04 


11.49±0.04 


10.63±0.05 


9.85±0.03 


7.86±0.02 



above process. The first row and last row are for MIR sources 
matched with suspected YSOs from Robitaille. The source 
in the middle row has no SIMBAD entry. The top source has 
the IRAC432 false colours of a PN but a star is seen at the 
same position in the 2MASS and optical images. The bot- 
tom two sources are in complex MIR regions close to bright 
MIR sources where diffraction spikes are playing a role in 
adversely affecting the GLIMPSE-I point-source photome- 
try. In Fig. [5] we give images for three of the 33 remaining 
high-quality PN candidates (excluding the two known and 
four newly confirmed PNe). The top two are associated with 
suspected YSOs in Robitaille et al. (2008) while the bottom 
source has no SIMBAD entry. 



After removing the two known PNe and four newly con- 
firmed PNe, the remaining 27 high-quality PN candidates 
(39% of the original sample) and their associated NIR and 
MIR photometry are presented in Table [XT] We believe the 
majority of these sources are likely to be PNe. 



5 SPECTROSCOPIC FOLLOW-UP 

Until we have further probed the capability of our MIR tech- 
nique to reliably reveal optically invisible PNe, we have to 
rely on follow-up optical spectra to validate MIR-selected 
candidates as PNe. Consequently, the four MIR-selected 
PN candidates with optical counterparts were observed, to- 
gether with PM 1-104, with the double-beam WiFeS image 
sheer integral field unit (Dopita et al. 2007) on the ANU 
2.3m telescope at Siding Spring Observatory in July 2011 
and June 2012. This powerful instrument provides simul- 
taneous red and blue arm spectra with a dichroic sending 
light to two detectors of similar design equipped with differ- 
ent gratings. WiFeS is an areal spectroscopic facility across 
a 36x25 arcsecond region of sky ideal for resolved sources. 
In the moderate seeing of SSO this instrument also ensures 
that all the flux from compact sources can still be acquired. 

Table [4] summarises the spectroscopic observations, 
gives the source IRAC designations from the GLIMPSE-I 
archive and the coordinates (accurate to 0.3 arcseconds). 
All four PN candidates have NIR and MIR diameters close 
to 3 arcseconds and are compact optically. 
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Legend: 

+ : Normal stars 
*: O-rich AG B stars 
O : Optical Carbon stars; 
A : Embedded Carbon stars; 
□ : Hyperluminous objects 
and OH/I R stars 
x: Exotica (TTau stars, 
reflection nebulae) 
• : UCHIIs. 



Note reddening vector 
in top left hand corner 



in 

i 

to 




-0,5 



[5,8)-[8,0] 



Figure 3. Detail of the IRAC [3.6] - [4.5] versus [5.8] - [8.0] colour plane with the positions of our four newly confirmed MIR selected 
PN over-plotted as red circles. The serendipitously confirmed MIR PN PM 1-104 is also plotted as a blue circle and falls outside our 
MIR selection box. The plot is adapted from Fig. 9 of Cohen et al. (2011). The general locus of points comprises the combined generic 
locations of 87 different types of IR point sources. Two rectangles are plotted for the median ±3 sem boxes for diffuse HII regions (dH) 
and for the colours of the entire known PN sample (P). 




* 



Figure 4. NIR 2MASS J,H,Ks, optical SHS and MIR IRAC432 false colour 1.5x1.5 arcmin images of a selection of three of the 37 
rejected PN candidates from the original 70. Images in the first and last rows are for MIR sources G313. 1931-00.2925 and G344.9571- 
00.0554 matched with YSO? stars in Robitaille et al. (2008). The middle-row source at G306. 6949-00.08213 has no SIMBAD entry. The 
top row shows a source with the IRAC432 false colours of a PN but the optical images reveal a star at the same position in 2MASS 
and the optical. The bottom two rows show candidates in complex MIR regions close to bright MIR sources where diffraction spikes are 
playing a role in confusing the GLIMPSE-I point-source photometry. 
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Figure 5. MIR 2MASS J,H,Ks, optical SHS and IRAC432 false colour 1.5x1.5 arcmin images of three of the 27 remaining high quality 
PN candidates that satisfy our environmental criteria based on examination of their 2MASS, optical and MIR images. Top: G308.4962- 
00.1623, mid: G305.3441-00.7489, bottom: G323.9540+00.4384. See Table [XT] for further details. 



Table 4. Summary of spectral observations performed with the MSSSO 2.3m and the WiFeS spectrograph. 



Object 


Date observed 


Dispersion 


Wavelength coverage 


Resolution 


Exposure 


Comment 






lines/mm 


Blue (A) 


Red (A) 


A pixel" 


- 1 B/R 


(seconds) 




GLIPN1823-1133 


03 July 2011 


1530B, 1210R 


4184-5580 


5294-7060 


0.36 ; 


0.45 


2x600 


Nod-and-Shuffle 


GLIPN1530-5557 


04 July 2011 


1530B, 1210R 


4184-5580 


5294-7060 


0.36 ; 


0.45 


1200 


Nod-and-Shuffle 


PM 1-104 


05 July 2011 


1530B, 1210R 


4184-5580 


5294-7060 


0.36 ; 


0.45 


600 


Nod-and-Shuffle 


GLIPN1557-5430 


05 July 2011 


1530B, 1210R 


4184-5580 


5294-7060 


0.36 ; 


0.45 


1800 


Nod-and-Shuffle 


GLIPN1642-4453 


18 June 2012 


708B, 1210R 


4184-5580 


5760-7030 


0.17 ; 


0.45 


2x2000 


Poor conditions 



5.1 The optical spectra 

The WiFeS data extraction was performed by one of us 
(MS) using the WiFeS data reduction pipeline (Dopita 
et al. 2010). Wavelength calibration was via standard arc 
lamps and spectrophotometric standard stars LTT9239 and 
HR8634 were repeatedly observed enabling flux calibration. 
The "nod-and-shume" observing mode was employed for 
the first three PN candidates. The telescope is nodded "off 
source" repeatedly to a nearby sky-region while the charge 
is shuffled to a different CCD area permitting excellent sky- 
subtraction though with an exposure time and noise penalty. 
Due to poor conditions during observations of the fourth 
candidate standard 'stare' exposures were taken. 

5.2 Results 

Fig. [6] presents three of the observed PN candidates with 
useful blue and red spectra. For GLIPN1823-1133 (top row 
Fig. 3) both blue [OIII] lines are visible but H/3 is not seen 
due to heavy extinction. The high [Nil] /Ha ratio in the red 



(-2.23) rules out a HII region (e.g. Kennicutt 2000). The 
observed [OIII] lines are at least a factor of 10 weaker than 
the [Nil] and Ha lines. The [SII] lines are well detected, 
permitting an electron density estimate of —2000 cm -3 . 

For GLIPN1530-5557 (mid row) only the brighter of the 
[OIII] lines is seen while in the red Ha is strong with only a 
small trace of [Nil] visible. This red spectrum is typical of 
high excitation PN though this is ruled out by the detection 
of Hel at 7065A. A PN ID is still strongly favoured. 

The bottom row gives spectra of GLIPN1557- 
5431(PM 1-104). Again only the stronger of the [OIII] lines 
is seen in the blue. The red spectrum is similar to that for 
GLIPN1823-1133 with a high [Nil] /Ha ratio of -2.11, also 
eliminating a HII region. The observed [SII] line ratio gives 
an electron density of —4200 cm -3 . 

Fig. [7] presents the WiFeS red spectra for GLIPN1557- 
5430 and GLI1642-4453 (no useful blue data were obtained). 
For GLIPN1557-5430 Ha and [Nil] are seen with the ratio 
0.59. This is at the high end of that found in HII regions 
(Kennicutt 2000) though the source is not diffuse and a PN 
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ID is indicated. The [SII] lines are also visible and provide 
an electron density estimate of ^2200 cm -3 . This is a very 
faint, compact source in the optical. For GLI1652-4453 the 
S/N is low but [Nil] and Ha are clearly detected in the ratio 
0.77 which eliminates confusion with a HII region, strongly 
confirming PN status. The [SII] again are visible and like- 
wise give an electron density estimate of ^1300 cm -3 . These 
emission-line spectra confirm the highly likely PN nature of 
all five observed candidates including PM 1-104. 

The observed line fluxes, ratios and derived elec- 
tron densities, n e , from our flux- calibrated WiFeS data 
are summarised in Table \E\ Fluxes are in units of 
10 -15 ergcm -2 s -1 A -1 and are not corrected for reddening. 
Entries with a trailing ":" indicate indicative values due to 
low S/N. None of the PNe have detectable H/3 emission to 
determine a Balmer decrement, but we can determine a use- 
ful upper limit to the H/3 flux based on the observed blue 
S/N and hence find a lower limit to the reddening. We also 
give an upper limit to the reddening for two PNe from an 
inferred lower limit to the H/3 flux, using our [O in] A5007 
flux and assuming a maximum A5007/H/3 ratio of ^20 (e.g. 
Acker et al. 1992). For two PNe we estimate reddenings from 
a comparison of the integrated Ha fluxes with available 6 cm 
radio fluxes (see Bojicic et al. 2011). 

Upper-limit E(B — V) estimates are also given from 
Schlafly & Finkbeiner (2011), who recently updated the 
Schlegel, Finkbeiner & Davis (1998) dust maps along each 
source sight-line. We refrain from quoting a final single ex- 
tinction estimate for each PN due to the different estimates 
used and their variation. This is the appropriate approach 
given the uncertainties though it is clear, as expected, that 
the extinctions are high. The quoted heliocentric line veloc- 
ities are based on application of the IRAF emsao package 
to the higher-resolution, red-arm WiFeS spectra, and are 
accurate to ~10kms -1 . We also quote integrated Ha fluxes 
measured from the SuperCOSMOS Ha Survey (Parker et al. 
2005) following the method of Gunawardhana et al. (2012). 
These Ha fluxes are amongst the faintest yet determined for 
any PN (see Frew et al. 2012). 

Details of the five spectroscopically observed sources 
are listed in Table [6] We adopt a new nomencla- 
ture for confirmed, MIR discovered, PN candidates as 
GLIPNhhmmiddmm of similar form to the MASH PN 
nomenclature (e.g. Parker et al. 2006). GLI indicates origin 
in the GLIMPSE-I footprint, PN indicates the object is a 
confirmed PN and the concatenation of the J2000 positions 
to hhmm and iddmm follows. Co-ordinates for PM 1-104 
have been updated from the incorrect published value. The 
PN status column reflects the same format used in MASH 
where T: true PN and L: likely PN. 



6 SUMMARY AND FUTURE WORK 

We have investigated the potential of the available MIR sur- 
vey data as a tool to uncover new PN candidates that would 
be hard or impossible to locate optically. The motivation is 
to develop robust MIR PN candidate selection techniques 
that can uncover the significant numbers of Galactic PNe 
hidden behind extensive curtains of dust. For this pilot study 
six MIR colour-colour selection criteria were applied to the 
GLIMPSE-I point source archive. These are based on the 



median values of the unique MIR colours of the 136 previ- 
ously known PNe that fall within the GLIMPSE-I footprint 
(and assumed representative of the overall Galactic PN pop- 
ulation). Only 70 candidate sources were returned. Most 
Galactic PNe are well resolved (e.g. only 5.5% of MASH 
PNe are compact/star- like) and so will not be found in the 
GLIMPSE-I point-source archive which also has a very re- 
stricted Galactic latitude coverage. These factors substan- 
tially reduce the number of obscured PN candidates found. 

Multi-wavelength image montages of each candidate 
were examined and four with faint optical detections in 
the SHS survey were found. Spectroscopy confirmed their 
likely PN nature. This result represents a clear validation of 
our general MIR selection technique to identify high quality 
PN candidates. This is because apart from their faint opti- 
cal signatures (due to the extinction not being too severe) 
they simply fulfil the MIR selection criteria we have devel- 
oped to identify PN candidates once their MIR images have 
been checked. We also confirm the PN nature of PM 1-104, 
serendipitously uncovered close to one of our MIR selected 
sources, and update erroneous positions for both PM 1-104 
and K 3-42 that fall in our sample. 

We demonstrate that false-colour images of MIR se- 
lected PN candidates are of high diagnostic value. They en- 
able the environmental context of the MIR point sources to 
be evaluated, showing that GLIMPSE-I point-source pho- 
tometry cannot always be taken at face value. We thus 
rejected 37 (54%) of the 70 MIR selected candidates as 
contaminants due to adverse MIR background effects, as- 
sociated dubious photometry or the de-blending of diffrac- 
tion spikes around bright stars into multiple spurious point- 
sources. In some cases the character of complementary 
mult i- wavelength optical and NIR data led to rejection. This 
left 27 high-quality PN candidates not including the four 
new PN confirmations and the two previously known PNe 
returned by the search. These results highlight the dangers 
of using GLIMPSE-I point source photometry in isolation. 

In future we will extend the MIR selection to the 
GLIMPSE-II and GLIMPSE-3D surveys and expand the 
search to resolved sources. Sky-background following, appli- 
cation of a threshold a given (low) percentage above this 
background and then running of pixel-connectivity algo- 
rithms will be used to isolate resolved but discrete MIR 
sources in a process directly analogous to that used on opti- 
cal data (e.g. Hambly et al. 2001). This should find resolved 
MIR sources and yield their integrated MIR magnitudes. 
Ultimately we plan to extend our MIR colour-colour tech- 
niques to the all-sky coverage of WISE which now enables 
alternative MIR false-colour images to be constructed. The 
WISE 3.4/xm and 4.6/im bands are directly equivalent to 
the first two IRAC bands at 3.6/xm and 4.5/xm. The final 
two IRAC bands at 5.8 and 8/xm do not have any direct 
WISE equivalent (with the closest being the WISE 12/rni 
band) though the WISE 22/mi band is similar to the MIPS 
24/im band. WISE can be used as a substitute for IRAC 
outside of the GLIMPSE regions with excellent sensitivity 
but poorer resolution (ranging from six arcseconds for the 
shorter wavelength bands out to ^12 arcseconds at 22/xm). 

If we can show that the sensitivity and resolution of the 
WISE MIR bands can provide the same diagnostic capability 
as for IRAC, then we can search for MIR PN candidates 
across the entire sky using essentially the same selection 
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Figure 6. Extracted blue and red 1-D spectra from the WiFeS data-cubes for GLIPN1823-1133, GLIPN1530-5557 and GLIPN1557-5431 
(PM 1-104). The PN spectral signatures are clear. Red and blue emission lines are identified in the top panel. See text for details 
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Figure 7. The final extracted red 1-D spectra from the WiFeS data-cube for confirmed PNe GLIPN1557-5430 and GLIPN1642-4453. 
There were no useful data obtained from the blue arm for these two candidates due to the heavy extinction. See text for further details. 



Table 5. Observed integrated line fluxes and line ratios as measured from our flux-calibrated spectra (summarised in Table [JJ. The 

1 and are not corrected for reddening (refer to the text for further details). 



fluxes are expressed in units of 10 15 erg cm 2 s 1 



Line 


Wavelength 


GLIPN1530-5557 


PM 1-104 


GLIPN1557-5430 


GLIPN1642-4453 


GLIPN1823-1133 




(A) 


F(A) 


F(A) 


F(A) 


F(A) 


F(A) 


H/3 


4861 


<0.3 


<0.5 


<0.2 


<0.1 


<0.4 


[O in] 


4959 


0.3: 


<0.5 


<0.2 


<0.1 


1.1: 


[O in] 


5007 


0.8: 


1.1: 


<0.2 


<0.1 


2.8 


He I 


5876 




0.5: 






0.5 


[O ii] 


6300 




2.0: 








[O ii] 


6363 




1.0: 








[N ii] 


6548 


0.3: 


15.6 


0.3: 


0.2: 


7.8 


Ha 


6563 


13.9 


31.5 


1.7 


1.16 


14.6 


[N ii] 


6584 


0.8 


51.0 


0.7 


0.67 


24.8 


He I 


6678 




0.4 






0.2 


[S ii] 


6717 




1.6 


0.1: 


0.13 


1.6 


[S ii] 


6731 




2.8 


0.15 


0.16 


2.3 


He I 


7065 


1.3 


1.8 




no coverage 


no coverage 


[N II] /Ha 




0.07 


2.11 


0.59 


0.77 


2.23 


[S II] 6717/6731 






0.58 


0.68 


0.80 


0.70 


n e (cm -3 ) 






4200 


2200 


1300 


2000 


E(B - V)n a /H/3 




2.8 - 4.5 


>2.9 


>1.0 


>1.3 


2.4 - 3.4 


E(B - V) 6cm/Hc , 






3.0 






2.7 


E(B - V)( dust _ map ) 




7.5 


4.2 


4.4 


5.2 


3.8 


A v 




23.3 


12.9 


13.6 


16.1 


11.7 


Vhel (kms -1 ) 




-21±10 


-63±10 


±20±10 


-52±10 


±105±10 


logF(Ha) (WiFeS) 




-13.86 ±0.10 


-13.50 ±0.05 


-14.79 ±0.15 


-14.94 ±0.15 


-13.84 ±0.10 


logF(Ha) (SHS) 




-14.2 ±0.2 


-13.46 ±0.10 


-14.9 ±0.2 


-15.3 ±0.2 


-14.2 ±0.2 



Table 6. Summary details of the new, confirmed MIR selected PNe. 



GLIMPSE-I ID 


PN-ID 


RA 


Dec 


I 


b 


PN 






(J2000) 


(J2000) 


degrees 


degrees 


status 


SSTGLMA G323.9366±00.2783 


GLIPN1530-5557 


15 30 41.6 


-55 57 27 


323.9365 


0.2784 


L 


SSTGLMA G327.8259-00.8710 


GLIPN1557-5430 


15 57 15.5 


-54 30 07 


327.8262 


-0.8712 


L 


SSTGLMA G327.8293-00.8879 


GLIPN1557-5431* 


15 57 21.0 


-54 30 46 


327.8292 


-0.8878 


T 


SSTGLMA G339. 7362-00.8468 


GLIPN1642-4453 


16 42 22.0 


-44 53 35 


339.7365 


0.8467 


T 


SSTGLMA G019.5325±00.7308 


GLIPN1823-1133 


18 23 59.9 


-11 33 39 


19.5326 


0.7309 


T 



* previously known as PM 1-104 
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criteria. Examination of known PN detected in WISE (e.g. 
Fig. |5| reveals potential in this regard. In this way we can 
compile MIR-selected PN candidates across the entire area 
covered by the SHS and IPHAS surveys and also to higher 
latitudes where there is no Ha coverage. This work is now 
underway. 
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